The density dependence of the transition temperature in a homogenous Bose fluid 
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Transition temperature data obtained as a function of particle density in the *He-Vycor system are 
compared with recent theoretical calculations for 3D Bose condensed systems. In the low density 
dilute Bose gas regime we find, in agreement with theory, a positive shift in the transition tempera- 
ture of the form AT/Tq — 7(na"^)^'^. At higher densities a maximum is found in the ratio of Tc/Tq 
for a value of the interaction parameter, na^, that is in agreement with path-integral Monte Carlo 
calculations. 
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The role of interparticle interactions in the determina- 
tion of the properties of low density Bose-Einstein con- 
densed (BEC) systems has been a topic of interest for 
many years. In spite of a long history of theoretical 
investigation dating back to the 1950s, elementary 
questions such as the possible shift in the transition tem- 
perature, Tc, with density and interaction strength have 
remained unsettled until the recent past. In the case of 
the repulsive interactions in the dilute Bose gas, there 
has now emerged a consensus Q [^ Q that Tc will be 
an increasing function of the interaction parameter, na'^, 
where a is the hard sphere diameter and n the particle 
density. This may seem a surprising result, since it is 
well known that in the case of liquid *He the superfluid 
transition occurs at a temperature well below the transi- 
tion temperature, Tq, for an ideal Bose gas with the same 
particle mass and number density. Moreover, a number 
of the earlier calculations [|| had found that the transi- 
tion temperature would be reduced as a consequence of 
interparticle interaction. 

Motivated by the recent theoretical work in this area, 
we have examined the dependence of the transition tem- 
perature on superfluid particle number for the *He-Vycor 
system. In our early work with this system we demon- 
strated, for the first time in 1983 g], an experimental re- 
alization of the weakly interacting or "dilute" Bose gas. 
The lowest density achieved in these experiments was on 
the order of 2 x 10^^ per cm^. This is sufficiently low to 
provide a region of overlap, in terms of the interaction 
parameter, with the values of na"^, currently accessible 
to the BEC systems realized with ^^Rb or ^^Na |] 
atoms confined within magnetic or optical traps. In the 
case of Bose condensed atomic hydrogen ^ the small s- 
wave scattering length and limits on the particle density, 
set by recombination, restrict this system to values of na'^ 
several orders of magnitude below the values that can be 
achieved with Bose-condensed Na, Rb or He. 

For questions such as the effect on Tc of increasing the 
interaction parameter, the *He-Vycor system offers ad- 
vantages over the BEC systems of trapped atomic gases, 
because in the Vycor case, the interaction parameter can 
be varied continuously from the low-density, weakly inter- 
acting limit to the strongly interacting regime currently 



inaccessible with the alkali vapor systems. Further, work- 
ing with the He- Vycor system allows much larger sample 
sizes, on the order a cm^, with an essentially unlimited 
time for observation of the BEC superfluid state. 

In this letter we will give only a brief summary of our 
experimental methods. The cryogenic techniques and 
thermometry methods are those standard in the study of 
superfluid ^He, and the reader is referred elsewhere px| ] 
for details. For the present experiments, we have used a 
torsional oscillator technique to obtain a signal propor- 
tional to the superfluid particle density within our Vycor 
sample, which then allows us to estimate the number den- 
sity for the particles participating in the superfluid and 
as well as to determine the superfluid transition temper- 
ature. The interior channels of the porous Vycor glass 
used for these measurements range in diameter from 4 
to 8 nm and form a highly interconnected 3D network. 
The superfluid helium atoms are constrained by van der 
Waals forces to move over the complex 3D-connected sur- 
face provided by the pores. It is important to appreciate 
that at low temperatures the thermal wavelength of the 
mobile superfluid particles can be larger than the pore 
size and that the Feynman exchange cycles characteriz- 
ing the BEC or superfluid state will link many pores at 
low particle density [Q. Therefore we model the su- 
perfluid phase as a homogeneous Bose gas constrained 
within the volume of the Vycor sample. We anticipate, 
however, that the influence of the substrate and pore ge- 
ometry may be reflected in a small modification of the 
effective mass of the mobile superfluid atoms. We then 
expect to observe, in the low density limit, thermody- 
namic properties similar to those of the "free" or ideal 
Bose gas with an effective mass, m*. 

In the discussion of our experimental results, we shall 
first consider the data for the dilute Bose gas regime. The 
quantities required for a comparison with theory are the 
particle number density and the corresponding transition 
temperature. The number density is estimated from an 
extrapolation of the superfluid signal to zero tempera- 
ture and the calibrated mass sensitivity of our torsional 
oscillator. Following the recent calulations M M H, we 
expect that for low densities T^, will be given by 



To[l + 7(na3)i/3] 



(1) 



where Tq = [2tiTt? /m*kBC,{i/2f/^]n'^/^ is the transition 
temperature for an ideal Bose gas with particle mass m* 
and density n. The coefficient for 7 is positive as required 
for Tc to rise above Tq as the interaction parameter, na^ 
increases. Although there is agreement on the form of 
Eq. 1, the theoretical estimates for 7 range over more 
than an order of magnitude, from an estimate of 0.34 
reported by Griiter, Ceperley, and Laloe (GCL) |^, 2.9 
by Baym et al. [Q, to a value of 4.66 given by Stoof [||. 
For a convenient comparison to experimental results 
we cast Eq. 1 in a linear form with the variable, n^'^a., 
as 



Tr 



277%^ 



n2/3 m*fcsC(3/2)2/3 



[l+7(ni/3a)]. 



(2) 



In Figure 1, we then plot Tc/n^/^ against the parameter 
n^/'^a, taking a value of 0.22 nm |l2| for the helium hard 
core diameter. As expected from theory, a linear fit gives 
a good representation for our data in the low density 
regime. The zero density intercept of the fit serves to 
determine an effective mass ratio of m*/m = 1.34, where 
m is the ^He mass, and the slope yields a value of 5.1 for 
7, which is within a factor of 1.09 of the value given by 
Stoof. 
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FIG. 1. The quantity Tc/n 
the parameter n^^'^ 



^" is plotted as a function of 
The straight line is our fit to these 



data. The arrows indicate the range of interaction parameter 
covered in the ^H, '*He, *^Rb, and ^^Na experiments. Note 
that in the case of ^H the scale has been enlarged by a factor 
of 10. 

It is of interest to compare the range of the interac- 
tion parameter explored in the various BEG experiments. 
The highest particle densities achieved to date are 4x10^^ 
atoms/cm'^ in the case of '^^Rb and 3x10^^ cm'^ for ^■^Na, 
however, for both the alkali atoms the s-wave scattering 
lengths are much larger than for the respective value for 
the ''He atom. The larger scattering lengths for the alka- 
lies more than compensate for the lower particle densities 



and result in a range of overlap in terms of the interaction 
parameter between all three systems. This is illustrated 
in Figure 1, where we indicate the range of the quan- 
tity n^/'^a for the '^He, experiments ^'^Na and ®^Rb ex- 
periments. Since the ^He-Vycor system clearly exhibits 
superfluid properties, one may also expect superfluidity 
to exists in the Bose-condensed alkali systems at compa- 
rable values of the interaction parameter and sufficiently 
low temperature. 
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FIG. 2. The observed transition temperature, Tc, is plot- 
ted as a function of the particle density, n. The top curve 
through the data is the function given by Eq. 1 with the pa- 
rameters for the effective mass, m* and 7 determined as from 
the linear fit shown in Fig. 1. The two other curves are for 7 
= 4.66 as given by Stoof and 7 = 0, the noninteracting case. 

In Figure 2 we show a more conventional plot of the 
low density data. Here we have plotted Tc as a function 
of the particle density, n. The curve through the data is 
the theoretical expression given in Eq. 1 with our best 
fit values for 7 and m* as determined previously. Two 
curves based on the 7 estimates of Stoof and Baym et al. 
are also shown. The lowest curve is for the case of the 
noninteracting Bose gas case (i.e., 7 = in Eq. 1). The 
relatively close agreement with the estimate of Stoof is 
clear in this plot. 

In contrast to the other theoretical treatments of the 
BEG problem for the interacting Bose gas H H], the 
calculations of GGL are not restricted to the low den- 
sity limit. The path- integral Monte Garlo technique em- 
ployed by GGL allows a prediction for the ratio of Tc/Tq 
over the entire range of na^, from the dilute gas limit 
to densities approaching the freezing point of liquid he- 
lium. In the top panel of Figure 3 we show the results of 
the GGL calculation taken from Ref. g. At low values 
of the interaction parameter, GGL also find that the ra- 
tio, Tc/Tq, increases in proportion to n^' '^a in agreement 
with the calculation of Stoof and of Baym et al. An in- 
teresting feature of the GGL calculation is the maximum 
in Tc/To found near a value of the interaction parameter 
na^ = 0.01. 

In the lower panel of Figure 3, we show the data ob- 



tained in a number of different ^He-Vycor experiments at 
Cornell, including a recent measurement designed specif- 
ically to map out the region of the peak found by GCL. 
Although the higher density data show more scatter than 
the data of Figs. 1 and 2, the trends are clear, and there 
is a pleasing agreement in the qualitative form of the 
experimental data as compared to the GCL calculation. 
In particular, there is an excellent match in the position 
of the peak value for the ratio Tc/Tq as a function of 
the interaction parameter. A remaining problem, how- 
ever, is the as yet unexplained disagreement in the mag- 
nitude of the effect estimated by GCL for low densities 
as compared to our experimental data or the calculations 
of Stoof § and Baym et al. |]. 
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FIG. 3. The top panel shows the values of the ratio Tc/Tq 
calculated by GCL as a function of the interaction parameter 
na . The curve through the lower density points is given by 
Tc/To = 1 + 0.34(na^)^''='. The curve through the higher den- 
sity data including the region of the peak is merely a guide 
for the eye. In the lower panel we plot the values of Tc/To, 
obtained from a number of different Vycor experiments (dis- 
tinguished by different symbols), as a function of na^. The 
dashed curve is the theoretical estimate of Stoof [2] . 

In conclusion, we wish to emphasis that studies with 
low density Bose condensed helium systems can provide 
a useful and illuminating contrast to the trapped gas sys- 
tems for the study of properties of the weakly interacting 
BEC fluids. Particularly important aspects of the ''^He- 
Vycor system are the ability to study long term phenom- 
ena, such as the question of robustness or stability of flow 
states in a Bose-condensed superfluid, and to extend ob- 
servations of BEC to the realm of larger interaction pa- 
rameters than are accessible to the trapped gas systems. 
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